Tachyon field in loop quantum cosmology: An example of traversable singularity 
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Loop quantum cosmology (LQC) predicts a nonsingular evolution of the universe through a bounce 
in the high energy region. But LQC has an ambiguity about the quantization scheme. Recently, 
the authors in [Phys. Rev. D 77, 124008 (2008)] proposed a new quantization scheme. Similar 
to others, this new quantization scheme also replaces the big bang singularity with the quantum 
bounce. More interestingly, it introduces a quantum singularity, which is traversable. We investigate 
this novel dynamics quantitatively with a tachyon scalar field, which gives us a concrete example. 
Our result shows that our universe can evolve through the quantum singularity regularly, which is 
different from the classical big bang singularity. So this singularity is only a week singularity. 

PACS numbers: 04.60.Pp, 04.60.Kz, 98.80.Qc 



I. INTRODUCTION 

Among the candidates for a theory of quantum gravity, 
the non-perturbative quantum gravity develops rapidly. 
In recent years, as a non-perturbative quantum gravity 
scheme, Loop Quantum Gravity (LQG) shows more and 
more strength on gravitational quantization [1]. LQG is 
rigorously constructed on the kinematical Hilbcrt space. 
Many spatial geometrical operators, such as the area, 
the volume and the length operator have also been con- 
structed on this kinematical Hilbert space. The success- 
ful examples of LOG include quantized area and volume 
operators 0, 0, [j, 0] , a calculation of the entropy of 
black holes Loop Quantum Cosmology (LQC) [7], etc. 
As a successful application of LQG to cosmology, LQC 
has an outstanding and interesting result — replacing the 
big bang^spacetime singularity of cosmology with a big 
bounce 8]. In addition, LQC also gives a quantum 
suppression of classical chaotic behavior near singular- 
ities in Bianchi-IX models 0, [l3|- Furthermore, it has 
been shown that the non-perturbative modification of the 
matter Hamiltonian leads to a generic phase of inflation 

Recently the authors in [l4| proposed a new quantiza- 
tion scheme (we will call it hms scheme in the following) 
for LQC [131 ■ In this new quantization scheme, the clas- 
sical big bang singularity is also replaced by a quantum 
bounce. The most interestingly, this new scheme intro- 
duces a novel quantum singularity. At this quantum sin- 
gularity the Hubble parameter diverges, but the universe 
can evolve through it regularly, which is different from 
the case for the classical big bang singularity. In order 
to investigate the novel quantum singularity in this (j-ms 
scheme, we apply it to the universe filled with a tachyon 
field and compared the result with the classical dynam- 
ics and the fl scheme which is presented in our previous 
paper [l3i] . 

The organization of this paper is as follows. In Sec. 



im we give out the effective framework of LQC coupled 
with the tachyon field, and a brief review of the new 
quantization scheme suggested in p^ . In addition we 
show some general properties of this specified effective 
LQC system. In Sec. |Tnl we use the numerical method 
to investigate the detailed dynamics of the universe filled 
with the tachyon field in the fiMS scheme. At the same 
time we compare the difference between the hms scheme, 
the p, scheme and the classical behavior. In Sec. IIVI we 
present some comments on the traversable singularity. 
In Sec|V]we conclude and discuss our results on the hms 
quantization scheme. Throughout the paper we adopt 
units with c = G = h = 1. 



II. THE EFFECTIVE FRAMEWORK OF LQC 
COUPLED WITH TACHYON FIELD 

The tachyon scalar field arises in string theor y [TBI , [l^ , 
which may provide an explanation for inflation jl7| at the 
early epochs and could contribute to some new form of 
the cosmological dark energy [l^ at late times. More- 
over, the tachyon scalar field can also be used to inter- 
pret the dark matter ^19,] . In this paper we investigate 
the tachyon field in the effective framework of LQC. Ac- 
cording to Sen [l^, in a spatially flat FRW cosmology 
the Hamiltonian for the tachyon field can be written as 



5H2 



where H^ 



is the conjugate momentum for the 
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tachyon field (f>, V{(f>) is the potential term for the tachyon 
field, and a is the FRW scale factor. So we have 

-1 < < 1. 

Following our previous work [isj . we take a specific po- 
tential for tachyon field [T6| in this paper, 

where Vq is a positive constant and a is the tachyon mass. 
Similar to [Q, we set Vq = 0.82 and a = 0.5. For the fiat 
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model of universe, the phase space of LQC is spanned by 
coordinates c = and p = a? , being the only remaining 
degrees of freedom after the symmetry reduction and the 
gauge fixing. In terms of the connection and triad, the 
classical Hamiltonian constraint is given by [l^l 



8717' 



(1) 



Considering ^ quantization scheme, the effective Hamil- 
tonian in LQC is given by |2lj 



H, 



■^sin^ (^c) + 



(2) 



The variable n corresponds to the dimensionless length 
of the edge of the elementary loop and is given by 



(3) 



where ^ is a constant > 0) and depends on the par- 
ticular scheme in the holonomy corrections. ^ is given 

by 



(4) 



where Ip is the Planck length. Considering the ^ms 
quantization scheme, the effective Hamiltonian in LQC 
is given by [iJl 



1 



H, 



1 



■ sin^(/xc) 



-Vp + H^- (5) 



In analogy to the jjq scheme and the p, scheme, we can 
understand this hms scheme as follows. We write the 
effective Hamiltonian as 



H, 



While fiMS is given by 



^sin^ (pMSc) + H^. (6) 



sm 



(PMSC) 



1- a/1 



■ sin (/^c) 



3/i 



(7) 



with p, determined by Eq. ([3]). As to the argument for 
the above equation we refer to [l3|- With notation of 
the effective energy density Pe// and the effective pres- 
sure Peff, we can write the modified Friedmann equa- 
tion, modified Raychaudhuri equation and the conserva- 
tion equation in the same form 



a 



Stt 



-Peff, 



(8) 

(9) 
(10) 



In the above equations, H = ^ stands for the Hubble 
parameter; pe// and Pg// are the effective energy density 
and the effective pressure, respectively. For the following 
three situations: the classical cosmology, the ft scheme 
and the fJ-MS scheme in LQC, we have 



Peff,cl 
Peff,cl 

Peff,fJ.MS 



V 



Pc J Pc 



p±_ 
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(11) 

(12) 
(13) 

(14) 
(15) 



I 1 - 4^ 



P^\ P4, 
3pc 



1 - ^ 
3 I 3pc 



3 , 
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3 Pc j i Pc 



(16) 



where pc = ■ In the first two lines we have al- definitions of the energy density and the pressure [l^l 

ready used the Hamiltonian for the tachyon field and the ^ ^tt 

p^^ — = . — (17) 
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Correspondingly, we have the following evolution equa- 
tions 



classical 



= - 1- 



T3.^ (l - 
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3 Pc 
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(18) 
(19) 

(20) 



In the above equations, "— " corresponds to the expand- 
ing universe while corresponds to the contracting 
universe. For the p, schenie, the bounce happens at 
Pc/, — Pc so we have [H, [2l|, \23^ 



< 



1 



For the um^ scheme, the bounce happens at — 3pc, 
so we have [l4| 



Voe-°"f' 



< 3pc 



The comparison of these features are shown in Fig. [T] 
and[2l In the left panel of Fig. [T]we compare the differ- 
ent behavior of the Hubble parameter versus the energy 
density of the matter. The upper half corresponds to 
the expansion stage of the universe, and the lower half 
corresponds to the contraction stage. For the fl quanti- 
zation scheme, we can see clearly the bounce behavior at 
P<t> = Pc- When p0 > pc/2 the universe meets a super- 
inflation phase {H > 0). When becomes small, the 
universe behaves as the standard picture. While for the 
Pais scheme, the bounce happens at p^ = 3pc- When 
> 3pc/2, the universe meets a superinflation phase. 
In the region where is small, the Hubble parameter 
of the pms scheme is smaller than the one of the stan- 
dard cosmology. In this sense the universe of the /ij\/s 
scheme expands slower than the universe of the standard 
cosmology. 

In Fig. [5] we compare the different behavior of in 
the expanding universe (taking "-" in Eqs. (fT8|) - (pO|) ). in 
the phase space, i.e., (p-cj) space. Firstly, from the upper 
three subfigures, we can see that the quantum correction 
in the p, scheme of LQC changes the amplitude of classical 
(j) very small, while in contrast, the quantum correction 
in the pms scheme changes this amplitude significantly. 



Besides, the quantum correction in the p scheme changes 
obviously the line shapes for iso-0, while the quantum 
correction in the pms scheme changes the line shapes 
negligibly. This difference results from the different loca- 
tions of the bounce region. In fact, the quantum effect is 
the strongest in the bounce region for different quantum 
corrections. The region shown in these three upper sub- 
figures is near the bounce region for the p scheme while 
some far away from the bounce region of the pms scheme. 
Secondly, in the lower three subfigures, we compare the 
quantum correction from pms of LQC with the classical 
behavior. We see that the behavior of near the singu- 
larity region is changed completely, (j) diverges when the 
state approaches the singularity line (the contour line of 

P4> = ^Pc) except for one point (j),(t>^ = (^^^ Iff ' 

where (f) — a. Yet near the bounce region, the behavior 
of (j) is similar to the one of the p, scheme in its corre- 
sponding bounce region. 



III. QUANTITATIVE ANALYSIS OF THE 
COSMOLOGICAL DYNAMICS COUPLED WITH 
TACHYON FIELD 

In the original paper [3] , the authors provided qualita- 
tive analysis of the dynamics for LQC in the pms scheme. 
Take the advantage of the specific model of the universe 
coupled with the tachyon field, we can investigate this 
dynamics for LQC quantitatively, and compare the dif- 
ference between the classical, the p, scheme and the pMS 
scheme dynamics. We solve the equations (ITO)) and 
(PO)) with the Rung-Kutta subroutine. The result is pre- 
sented in Fig. [21 

For the pMS scheme, the difference between the classi- 
cal and the quantum behaviors is more explicit in the re- 
gion between the bouncing boundary and the singularity 
line. While the difference in the region on the right-hand 
side of singularity is negligible. Note that the admissible 
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states for the scheme all locate in this region. So we can 
imagine that the difference between the Hms scheme and 
the fL scheme is nothing but the difference between the 
classical one and the p, one, which is presented in Figs. 1 
and 2 of our previous paper [T3| . Near the bounce region, 
the quantum behavior is similar for both the /2 scheme 
and the fiMS scheme. Certainly, this behavior emerges 
at different places in space for the p, scheme and the 
Mms scheme respectively. Near the singularity region, 
there is not special respects for the quantum evolution 
in the space, except that the quantum trajectories 
are much steeper than the classical ones. This steeper 
behavior is resulted from the singularity behavior of the 
quantum dynamics which makes (p much larger than the 
original classical ones. From the right panel of FiglJl we 
can see that both the hyper-inflationary and deflation- 
ary phases of the universe emerged clearly. The universe 
expands increasingly faster before singularity untill the 
acceleration becomes infinity. This stage corresponds to 
the hyper-inflationary phase for the universe. After this 
singularity, the universe expands more and more slowly. 
The behavior comes back to the classical one [U, 
quickly. This stage is the deflationary phase for the uni- 
verse. 



IV. COMMENTS ON THE TRAVERSABLE 
SINGULARITY 

As a dynamical system, ((20)) is singular when = 
Spcfi, which corresponds to the right "C" shaped curve 
of the left panel of Fig. [3] That means that the dynam- 
ical system is only defined in two separated regions. One 
is the region between the left "C" shaped line and the 
right "C" shaped line, and the another is the region on 
the right-hand side of the right "C" shaped line. Then 
a question arises naturally-does the numerical behavior 
traversing the separated "C" shaped line make sense or 
is it only a numerical cheating [35]? If we consider these 
two regions separately the dynamical system is well de- 
fined in the sense of the Cauchy uniqueness theorem. The 
orbits in the phase space (j)-4> are smooth. Taking the 
phase space as R^, these orbits are smooth curves. So 
these smooth curves have proper limit points on the sep- 
arated "C" shaped line. From the numerical solutions 
shown in Fig. [31 we can see that the different orbits 
have different limit points on the "C" shaped line. This 
is because the vector flow generating the trajectories (on 
the 4>-(p plane) has well-defined directions (vertical) at 
the singularity, which means the trajectories cannot in- 
tersect there. Then it is natural to join the two orbits in 
the two separated regions with the same limit points. In 
this way we get a well-defined dynamical system in the 
total region which lies on the right-hand side of the left 
"C" shaped line. We can expect the numerical solution 
with the Rung-Kutta method will converge to this solu- 
tion (we also performed numerical integrations starting 
from both sides towards the singular point and obtained 



the same result). So the numerical result presented in 
the above sections does make sense. 

In the following, we come back to the spacetime to 
check the property of this traversable singularity. Cor- 
responding to the two separated regions of phase space 
(j)-(f), the scale factor a can be determined by (p a-nd 
through 

a(0,0) = (p^/i/^)i/^ (21) 

which is a smooth function of (j) and (f>. Therefore, the 
two spacetime regions of the universe corresponding to 
these two regions is smooth. When the universe evolves 
to — 3pc/2, a is also well defined through (PT|) . That 
implies that the whole spacetime of the universe is well 
defined through joining the above mentioned two regions 
together by the well-defined a at p^ = 3pc/2. But since 
(j) is singular there, a is also singular there, which makes 
the spacetime of the universe unsmooth. In this sense the 
spatial slice corresponding to this special universe time 
is a traversable sin gularity of the spacetime. 

According to [26l.l27l.l28l.l29t. singularities can be classi- 
fied into strong and weak types. A singularity is strong if 
the tidal forces cause complete destruction of the objects 
irrespective of their physical characteristics, whereas a 
singularity is considered weak if the tidal forces are not 
strong enough to forbid the passage of objects. In this 
classification the singularity discussed here is only a weak 
singularity. As to the cosmological singularities, they 
can be classified in more details with the triplet of vari- 
ables (a, p,P) [aS IM [32, [13: Big bang and Big Crunch 
(a = 0, p and curvature invariants diverge at a finite 
proper time); Big Rip or type I singularity (a, p, P and 
curvature invariants diverge at a finite proper time); Sud- 
den or type II singularity (a and p are finite while P di- 
verges); type III singularity (a is finite while p and P 
diverge); type IV singularity (a, p, P and the curvature 
invariants are finite while the curvature derivatives di- 
verge). For the singularity discussed in this paper, a, p 
and P are all finite because of the regularity of cj) and 4> 
at the singularity. But the Ricci curvature invariant, 

R = 6{H^ + -), (22) 
a 

diverges at the singularity. So the singularity here does 
not fall in any type of above clarification. But it is more 
similar to type IV singularity than to the other types. 

V. DISCUSSION AND CONCLUSION 

In the classical cosmology our universe has a big bang 
singularity. All physical laws break down there. LQC 
replaces this singularity with a quantum bounce and 
the universe can evolve through the bounce point reg- 
ularly. Considering the quantum ambiguity of quanti- 
zation scheme in LQC, the authors in [14| proposed a 
new scheme. In addition to the quantum bounce, a novel 
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quantum singularity emerges in this new scheme. The 
quantum singularity is different from the big bang singu- 
larity and is traversable, although the Hubble parameter 
diverges at this singularity. In this paper, we follow our 
previous work |13j to investigate this novel dynamics with 
the tachyon scalar field in the framework of the effective 
LQC. We analyze the evolution of the tachyon field with 
an exponential potential in the context of LQC, and obvi- 
ously, any other choice of potential can lead to the similar 
result. 

In the high energy region (approaching the critical den- 
sity pc), LQC in the new quantization scheme greatly 
modifies the classical FRW cosmology and predicts a 
nonsingular bounce at density 3pc, which is located at 
a different density region compared with the p, quantiza- 
tion scheme. In addition, besides this quantum bounce. 



this new quantization scheme also introduces a quantum 
singularity, which emerges at density 1.5pc- At this quan- 
tum singularity the Hubble parameter diverges. But this 
singularity is different from the classical one. The uni- 
verse can evolve through this quantum singularity regu- 
larly. Different from the [1 scheme, the dynamics of the 
new quantization scheme will deviates from the classical 
one even in a small energy density region (see Fig. 1). 
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FIG. 1: (color online) The dotted line, dashed line and solid line correspond to the classical, the fi scheme and the fiMS 
scheme dynamics respectively. Left panel: Hubble parameter versus the energy density of matter. When pif, <C 1, the quantum 
correction becomes negligible in the p, scheme, i.e., the behavior of universe goes back to classical one. But for the hms the 
quantum correction is strong even when p0 <C 1. Right panel: The phase portrait for all admissible and (p and the graph 
shows the region boundary. These boundaries are nothing but the contour lines of in — space. The dashed line, dotted 
line and solid line correspond to oo, 3pc and pc respectively. The admissible states for classical system are between the two 
dashed lines. The admissible states for quantum system are the right part of the corresponding lines. 
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FIG. 2: (color online) Comparison of <?!> in the expanding universe. The horizontal axis is 0, and vertical axis is (j). The lines 
are contour lines for 0. The upper row is comparison of three dynamical systems near the bounce region of the p quantization 
scheme. The lower row is comparison of the classical system and the pji/s scheme. In the three subfigures of upper row we only 
show the admissible region for the p scheme, which is the common region for three dynamical systems. The lines are contour 
plot of (j). The figure does not show this singularity properly due to the numerical resolution. But the envelopes of islands in 
this figure are true. At present, we are not clear what physical phenomena these islands will result in. Near the bounce region, 
the behavior of is similar to the one of the p scheme in its corresponding bounce region (compared middle subfigure in the 
upper row with the one of last subfigure in the lower row) . 
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FIG. 3: (color online) Lept panel, comparing of the dynamical trajectories for the classical (dotted lines) and the /iMS quan- 
tization scheme (solid lines) under the expanding universe scenario in <j> — (p space. The left "C" shape line is the boundary 
of admissible states for the fiMs scheme, which is also the iso line of = 3pc- The right "C" shape line corresponds to the 
singularity (the iso line of — l.bpc) for the pMs scheme, which is also the interface for hyper-inflationary and deflationary 
states. Right panel, the concrete example of the dynamical behavior of and the Hubble parameter H near singularity for 
the fiMS quantization scheme, which corresponds to any line which passes through the singularity for the pMS scheme in the 
left panel. The behavior of other ones are similar on the singularity. The solid line shows behavior of 0, while the dotted line 
shows behavior of the Hubble parameter H = ^. This subflgure is a blow-up near singularity. 



